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1 Introduction 

Let {B^ 1 = (Bj, . . . , Bf),t > 0} be a <i-dimensional fractional Brownian motion (fBm) with Hurst 
index H in (0, 1). Let B 11,1 and B H ' 2 be two independent copies of B H . If Hd < 2, then the 
intersection local time of B H}1 and B H ' 2 exists (see [4]) and can be defined as 



a{h,t 2 )= I ' j 2 5{B^ 1 -B^ 2 )dudv, 
Jo J 



where 5 is the Dirac delta function. For any t\ and t 2 in R + , define 

X(t 1 ,t 2 ) = B»> 1 -B% 2 . 

We see that X = {X(ti,t 2 ),ti,t 2 £ K + } is a (2, (i)-Gaussian random field and satisfies the following 
scaling property: for any c > 0, 

{x(ct l ,ct 2 ),t l ,t 2 el + } = {c H x{t 1 ,t 2 ),t 1 ,t 2 eR+}. (1.1) 

If Hd < 2, then, for any x in M. d and rectangle E = [a±, b\] x [02, b 2 ] in the local time L(x, E) of 
X exists and is continuous in x, see [6]. When E = [0,t\] x [0,t2]> o>(ti,t 2 ) = L(0,E). Throughout 
this paper, we assume Hd < 2 to ensure the existence and the continuity of L(x,E). 

For any integrable function / : M rf — > K, one can easily show the following convergence in law 
in the space C([0, oo) 2 ), as n tends to infinity, 



n Hd-2 



/ / /(B^ 1 -B^ 2 )dudv, h,t 2 >0) A \a(h,t 2 ) / f(x)dx, t u t 2 >0\. 
J0 JO J L JR d > 
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In fact, letting E = [0, t{\ x [0, £2]) using the scaling property of the process X(u,v) in (1.1) and 
then applying the continuity of L(x,E), we get 



Hd-2 

n 



/ / f(B^-B^)dudv k n Hd / f(n H (B^-B^))dudv 
Jo Jo Jo Jo 

= n Hd [ f(n H x)L(x,E)dx 

= [ f(x)L(^,E)dx 
jV n H 



a(ti,t 2 ) / f(x)dx, 



where — denotes the almost sure convergence. 

If we assume J^ d f(x) dx = 0, then the random variable 

rnti rnt2 



n Hd-2 



/ / f(B^-B^ 2 )dudv 
Jo Jo 



converges in law to as n tends to infinity. It is natural to ask if there is a /3 > Hd — 2 such that 



/ / /(B^-B^dudv 
Jo Jo 



converges to a nontrivial random variable. This will be proved to be true. In order to formulate 
this result we introduce the following space of functions. Fix a number j3 £ (0, 2), define 

Hg = {/ G L l (R d ) : / |/(x)||a;|^dx<oo and / f(x)dx = o\. 

JR d jR d > 

For any / € Hq , by Lemma 4.1 in [3], the quantity 

11/11? = - / f(x)f(y)\y-x\-?dxdy 
is finite and nonnegative. The next theorem is the main result of this paper. 

— — d 

Theorem 1.1 Suppose < H < | and f £ Hq . Then, for any t\ and t% > 0, 

"nil rnt2 



n 2 



as n — > oo, where 



- / / f{B^-B^)dudv A v^dll/IU-dV^t^yC, 

f 27T") 2 Jo / 



(2-7r)2 Jo jo 

and C is a standard normal random variable independent of the processes B H ^ and B H ' 2 . 
In [3], Hu, Nualart and I proved the following functional central limit theorem 

[n^ £ f(B H (s))ds, t>0} -A {^C^ d \\f\\^_ d W(L t (0)),t>0}, 

2 



where W is a real- valued standard Brownian motion independent of B H and L t (x) is the local time 
of B . This paper can be viewed as an extension of the result in [3]. To prove our main result 
Theorem 1.1, we use the method of moments. Some techniques in [3] will be used, but new ideas 
are needed. The basic idea of the approach used in this paper is to apply the method of moments 
to a functional. When dealing with an integral on [0, ti] 2m x [0,t2] 2m 5 with respect to the measure 
du\ ■ ■ ■ du2mdv\ ■ ■ ■ dv2m , we make the change of variables W2k-i = n{u2k — "U2fc-i)i w 2k = U2k, 
S2k~i = n(v2k — V2k-i), S2k = V2k, I < k < m. Then, the increments of B H)l — B H ' 2 in small 
rectangles will be responsible for the independent noise appearing in the limit. This methodology 
could be applied to other examples of functionals and multi-parameter processes. 

Note that the constant Djj^ is finite for any H > We conjecture that our result is also 

true for < H < |, but we have not been able to show our result in the case H < The 
main reason is that we need to use Fourier analysis in the proof of our result. For example, we 
need to assume Hd > 2 — H in Lemma 4.2. 

In the Brownian motion case (H = ^ and d = 3), the functional version of Theorem 1.1 can 
be proved using a theorem by Weinryb and Yor [5]. A second order result for two independent 
Brownian motions in the critical case d = 4 and H = ^ was proved by Le Gall [2]. However, not 
nearly as much has been done for the case H ^ ^ and Hd = 2. The general asymptotic results for 
additive functionals of k independent Brownian motions were obtained by Biane [1]. This paper 
extends some results in [1] to fractional Brownian motions. General extensions are still largely 
unknown. 

After some preliminaries in Section 2, Section 3 is devoted to the proof of Theorem 1.1, based 
on the method of moments. Throughout this paper, if not mentioned otherwise, the letter c, with 
or without a subscript, denotes a generic positive finite constant whose exact value is independent 
of n and may change from line to line. We use l to denote 



2 Preliminaries 



Let {B^ 1 = (Bf, . . . ,Bf),t > 0} be a d-dimensional fractional Brownian motion with Hurst index 
H in (0,1), defined on some probability space (Q,,J-,P). That is, the components of B H are 
independent centered Gaussian processes with covariance function 

E(BiBi) = l -(t 2H + s 2H -\t-s\ 2H ). 

We shall use the following property of the fractional Brownian motion B H . 

Lemma 2.1 Given n > 1, there exist two constants c\ and C2 depending only onn, H and d, such 
that for any = sq < si < ■ ■ ■ < s n and X{ E M. d , 1 < i < n, we have 

n n n 

Cl \xi\ 2 {8i - Sl _i) 2H < Var(£> • {Bl - < c 2 ]T \ Xi \ 2 {s t - s^) 211 . 

i=l i=l i=l 

Proof. The second inequality is obvious. So it suffices to show the first one, which follows from 
the local nondeterminism property of the fractional Brownian motion; see, e.g., [1] and [3]. ■ 

The inequalities in Lemma 2.1 can be rewritten as 

n n 2 n n n 2 

^~ Si ^ 2H ^ Var (l>*-^) ^ c 2^|I>j| (si-s^) 211 . (2.1) 

i=l j=i i=l i=l j=i 
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The next lemma gives a formula for the moments of the random variable W a(ti,t 2 ) C appearing 
in Theorem 1.1. 



Lemma 2.2 For any pGN, 



E[v / «(ii>*2)C] ! 



(2m md" / » ( det A(u, v)) 2 dudv ifp = 2m, 



_ otherwise, 
where E = [0, t\] x [0, £2] and A(u, v) is the covariance matrix of the Gaussian random field 

(B^-B^,l<i<m). 

Proof. This follows easily from the properties the normal distribution and the intersection local 
time a(t±, £2). ■ 



3 Proof of Theorem 1.1 



By the scaling property of X(ti,t 2 ) in (1.1), we see that, as random variables, 



Hd-2 



nti pnt'2 



I I f{B^ 1 - B^' 2 )dudv = n 
/o Jo 



£ 2 + Hd 
2 



ti r-t 2 



f{n H (B^-B^))dudv. 



10 Jo 

Therefore, it suffices to show Theorem 1.1 for the random variable 

2+Hd f' 1 f' 2 

Fn (t 1 ,t 2 )=n— / f{n H {B^ - B^ 2 )) dudv. 
Jo Jo 

The proof of Theorem 1.1 will be done in two steps. We first show tightness and then establish 
the convergence of moments. 



3.1 Tightness 

Tightness will be deduced from the following result. 

Proposition 3.1 For any integer m > 1, there exists a positive constant C independent of n such 
that 

1 2m . „ r f , .2 



E[F n (t u t 2 )Y m <C / \f{x)f(y)\\y\«- d dxdy 



Proof. Note that 



E[F n (ti,t 2 ] 



1 2m. 



n 



2m 



m(2 +Hd ) E / y^ {B ^-B^))dudv 



(3.1) 



i=\ 



where E = [0,*i] x [0,t 2 ]. 
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Using Fourier analysis and making proper change of variables, 

„ 2m 

(2vrn") 2 ^E [ / \{f{n H {B^ - << 2 )) dudv 



E 2 



8=1 



imd J J7^2m 



2 fYl, 2 Til- 2 7T1 

fl /(*) exp { - ivar ( £ ft ■ (sj 1 - << 2 )) - ^ ducfo 

i=l i=l i=l 



i=l 

2/11 



/ / n/^iK^- 1 ; 

j 2m ^ 2m 

x exp { - - Var ( E ft ■ <) - 5 Var ( E ft ■ B^j } <fa cfo d£ tfe, 



(3.2) 



i=i i=i 

where in the last equality we used the fact that f^ d f(x) dx = 0. 

Let t = maxjii, £2} and 3? be the set consisting of all permutations of {1, 2, . . . , 2m}. Set 

„ -. 2m 

I t (0= / exp{--Var(^e r <; 



(in. 



For any cr £ define 



2 m 



if(0= / exp{-ivar(X>-<)} 



where D ff = {i*G [0,t] 2m : u ct(1) < • • • < u a(2m )}. 
Therefore, it(£) can be decomposed as 

^t(o = E w 

For simplicity of notation, set 

2m 



2m 



i=l 



i=l 



From (3.1), (3.2) and (3.3), we can write 

E[F n (t u t 2 )] 2m < Cl [ <S> n (t,z){l t (0) 2 dadz<c 2 Y] [ $n(t,z)(l?(0) 2 dt 
Observe that 

„ s 2m ^ 2m v 

{lt(0) 2 = / exp --Var(E^-<)-2 Vaj (E^-<)f^^ 

„ f 1 2 m 2m. n 

= / exp --Var(Eer<)-2 Vax (E^ >5 ?)r U;ds ' 
where D = \ w, s G [0, t] 2m : ioi < ■ ■ ■ < W2 m and si < • • • < S2m}- 



(3.3) 



(3.4) 



cfe. (3.5) 



Using the second inequality in (2.1), we obtain that 

/ (if (o) 2 de<fe 

J R 4md 

is less than or equal to 

2m 2m 

/ / $ n (£, z) exp { - ?f V | V e,f [K - w^f" + (* - Si _!) 2 ^] ) dw ds d£ 



dz, 



with the convention wq = sq = 0. 



i=l j=j 



2/u 



Making the change of variables ?jj = ^ £j, Ui = Wi — Wi-i and Vi = Si — Sj_i for i = 1, . . . , 2m 
gives 



2m 2m 



< / / in ri i ex p (^4 • - »»)) - ii 

Jv*»* J[o,t] im f=i f = V n 

2m 

x exp | - y^|i?i| 2 (jiP + uf)}(ij?^*(fz, (3.6) 



with the convention r/2m+i = 0. 

Let y/KjjXi, . . . , yJ^nXim be independent copies of the (f-dimensional standard normal random 
vector and X2 m +i = 0. Then inequality (3.6) can be written as 



< c 3 n m{2 ~ Hd ^ 



2 m 



2m 



n i/(^)i iik 2h + 



K 2md J[ 0;t ]4m . = 1 ^ 



2m 

n 



x ll| ex P [ L 7Jt 

i=l 



X; 



du (if dz 



(3.7) 



To make use of the independence of X±,X2, ■ ■ ■ , Xi m , we replace the terms 



Xi 



i+1 



Xi 



exp i— w ■ 

K n H - I 2H , ,,2fl 
"i+l + v i+l 



U 2H + 



i = 2,4, ... ,2m 



6 



on the right hand side of inequality (3.7) with 2 and then obtain 



* n (t,z)(i?(s)yd£dz 



m 

x | exp 



2m 2m 

n i/(^)i u( u i H + 

R2md J[ Q j]4m ^ 



1 = 1 



k=l 



X 2 k 



X 2 k- 



n 



H 



C4 n 



x E 



,2-Hd 



' U 2k+ V 2k 
2 2 



2H i „,2H L,2H , „,2H 



U 2k~l + V 2k-1 



du dv dz 



^J[0,rt=! i=i 

X 2 
= 

n 



uf 1 + J u\ H + up 



d«i dv,2 dv\ dv2 dz\ dz2 



< c 5 t 



2-Hd 



\f(z 1 )f(z 2 )\z 1 \H d dz x dz-< 



(3. 



where in the last inequality we used Lemmas 4.1 and 4.2. 
Combining (3.8) and (3.5) gives the desired inequality. 

3.2 Convergence of odd moments 

Let t = max{*i, t 2 }. For any p £ N, y £ R d and £ = . . . , £ p ) E (R d ) p , define 

P c 
p(2-Hd) -,—r . . y j jj . 

= n |/(%)(e-^ - 1)| 

and 



(3.9) 



= / / <M£: ?/) { - J Var ( V ■ - <' 2 ) ) ) d« dv d£ dy. 

Note that for p = 2m, $n,p(£jy) is precisely the function defined (3.4). Also in the proof of 
Proposition 3.1, we have that E [F n (ti, t 2 )] 2m < cH n ^ 2rn . We are going to use see that if p is odd, 
then H n , tP converges to zero as n tends to infinity. This will imply the convergence of odd moments. 



Proposition 3.2 If p is odd, then 



lim H n „ = 0. 

n— >oo ' 



Proof. Using similar notation as in the proof of Proposition 3.1, we have 

r r /• P P 

p(2-Hd) 

H n p < ci n 2 E 



n i/(w)i U( u i H+v " H r 

RP d J[0,t] 2 P fJl 



i=l 



n l ex p { l ^h 



i=l 



X 



i+1 



Xi 



u 2 ^ + v™ 



U 2H + y 2H 



du dv dy 
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with the convention X p+ i = 0. Since Xi,X 2 , ■ ■ ■ , X p are i.i.d, 



p(2-Hd) 

H n . P <c 2 n 2 E 



P+l — OJMUCAl,A2,'" ,A p 



2 * 1 t,F)"5 



p+1 

2 



exp 



fc=i 

(2-Hd) 

< C3 n 2 



y2fc-i 

H 



X 2 k 



X 2 k- 



n 



U 2k~l + V 2k-1 



ds dt dy 



2 d . . N~ 



l/(yi)/(j/2)M H d yi dy 2 



x E 



R d J[0,t] 2 



exp 



;? 



citip cfc p dy p 



(Hd-2) 

< C4 n 2 



|/(yi)/(2/2)|yi| H d dyidy 2 



\f(y P )\\y P \ H d dy P 



where in the last two inequalities we used Lemmas 4.1 and 4.2. Therefore, lim H np = 0. ■ 

Propositions 3.1 and 3.2 show that -£f„ )P is uniformly bounded in n. Moreover, Proposition 3.2 
implies the following convergence of odd moments. 

Proposition 3.3 Suppose p is odd, then 

lim E[F„(ti,t 2 )] p = 0. 

n— >oo 

Proof. Since \M(F n (ti,t 2 )) p \ < -£f n . p for all p, this follows immediately from Proposition 3.2. ■ 



3.3 Some technical lemmas 

To prove the convergence of even moments, we need some technical lemmas. 
Recall that 

2m 

1 2m 

L JE 2m 



E [F n (t 1> t 2 )] am = n^'+^E / ^ J] f( nH K A - nH K 2 ) dudv 

3=1 



where E = [0,<i] x [0,t 2 ]- 

Let ^ be the set consisting of all permutations of I = {1, 2, . . . , 2m} and 



jo = u < ui < u 2 < ■ ■ ■ < u 2m < ti, = v < v± < v 2 < ■ ■ ■ < v 2m < t 2 j. (3.10) 



Then 



2m 

E[F n (t u t 2 )] 2m = (2 m )\n^ H ^J2^[ Ylf{n H B^-n H B»? {j) )dudv\ (3.11 



For any e > 0, define 
H 



'n,2m,£ 



/ / $n, 2m (£, y) exp { - - Var ( £ & ■ - } d« dv <% dy 

J R 4md J Re v I j = l 



where 

Re = |0 < U\ < U2 < ■ ■ ■ < U2m < t, < V\ < V 2 < ■ • • < V 2m < tj 

P| ||«2j - Uj-2\ < e or \v2j - V2j-i\ < e> for some j £ {1,2, ... , m}| 
with the convention uq = v o = and t = maxjii,^}- 



Lemma 3.4 For any a G & , 



lirnsupH£ = 0. 



Proof. Note that 

Re = ^T=\ (Re,£,l U Re/,2) , 

where R e ,t,i = Re H {\u2£ — 1*20-2 1 < and R e ,e,2 = i? e H {|v 2 ^ — ^2^-2! < e }- So it suffices to show 
that 

lim sup / / <& n , 2m (£, y) exp { - - Var (V( r (fi^ 1 - ))} du dv d£ dy = 

for all I = 1, 2, ... , m and for i = 1,2. We will consider only the case i = 2 and the case i = 1 
could be treated in the same way. 

Let 

„ - 2m 

■MO = exp { - -Var ( £ & • Z^) } du 



J'=l 



and 

„ 1 2m 

■tf' e (0 = / exp { - -Var ( £ £ r < ) } 

•/{0<ltl<W2<-<W2m<t,|WM— "2^-2|<e} Z j = l 

Applying Cauchy- Schwartz inequality, we obtain 

r r 1 2m 

/ / $n,2m (£, y) exp { - -Var (V ■ (fl J 1 - ) ) } du dv d£ dy 
r 1 r 1 

< ( / ^nMC,y)(Jt(0) 2 ^d y y( / *„ |2 mfoy)(#' e (£)) 2 <*e<fo) 3 

<(#n,2m)H / ^n,2m(e,y)(«/r(e)) 2 ^^) 5 - 

Note that 

/ <S>n,2 m (ty)(Jr(0) 2 dtdy 

= / / $n, 2m (^ y) exp { - J Var (f? r (i^' 1 - i^' 2 )) } du dv d£ dy, 

J R 4md J Rc l 12 V ^ J 

where 

Re 1 = H Rel.2- 
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Since H n ^ m is uniformly bounded in n, we only need to show that 
r r i ^m 

lim / / $ n , 2m (e,y)exp{--Var(Ve, - - <' 2 ))) dn^d^y = 

for all £ = 1,2, ... ,m. 

Using Lemma 2.1 and then making the change of variables Wi = U{ — si = vi — Vi-\ and 

^2r 

2m 



rji = Yl^=i £j f° r * = 1)2,..., 2m with the convention uq = vq = and r\2m+\ = 0, we obtain 



$n,2m(£, y) exp { - ivar ■ (B*; 1 - <' 2 )) } dudv d$ dy 

„ „ 2m 2m 

< / / n n i ex p ('4 • - *)) - !i 

./r*™« J Re e AX XX n H 

2m 

x exp | - J^|?7i| 2 (u;P + s 2H )}dw (is dr/dy, (3.12) 



where 



t=i 



2m 2m 



.R^ = [0, t] 4m n I ^ < t, ^ Sj < t, i02< + W2i-i < e, S2^ + s 2 t-\ < e|. 



i=l i=l 

Using the same argument as in the proof of Proposition 3.1, we can prove that the right hand side 
of the inequality (3.12) is less than a constant multiple of e 2 ~ Hd . Letting e — > completes the 
proof. ■ 

For 1 < k < m, we define 

2m ,k = \u, v € [0, t] 2m :ui<u 2 <---< u 2m , v i < v 2 < ■ ■ ■ < v 2m , 

U2k ~ U 2k - 2 ^ V 2k ~ V 2k -2 ^ 1 

< u 2k - u 2k _i or < v 2k - v 2 k-ij- 

Recall the definition of & n ,2m(£,,y) in (3.9). The following result states that the integral over the 
domain 2m ^ k does not contribute to the limit of the 2m-th moment, which will play a fundamental 
role in computing the limits of even moments. 

Lemma 3.5 For any 1 < k < m, 

lim / / $ n , 2m (C,y)exp{--Var(^Ci-(<' 1 -<' 2 )))^^^(iy = 0. 



Proof. Define 



2m ,k = {u, v G [0, t] 2m :ui<u 2 <---< u 2m , vi <v 2 < ■■■ < v 2m , 

u 2k - u 2k _ 2 v 2k - v 2k _ 2 -, 
^ < u 2k - u 2k -i, < v 2k - V 2 k-\). 
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Using Cauchy-Schwartz inequality, we obtain 
r r i 2m 

/ / $n,2m(£, y) exp { - -Var (V & • ^ ~ K'^) } # dn du ^ 

JR 4md J 02m,k j = l 

/ 4 * n, 2m (e, y) / exp { - -Var ( £ & ■ (fl*' 1 - B% 2 )) } d w ^ dy ) 



< C! ( / $ n , 2m (C, y) / exp { - -Var ( £ £j ■ (B^ 1 - B^ 2 )) } d U d W d£ dy 

V jR 4md J 02m, k j = \ 



So it suffices to show 

2m 



lim / $ n ,2m(e, y) / exp { - Var ( V ■ (^J 1 - B*' 2 )) } du dv d V dy = 0. 

For j = 1,2,..., 2m, we make the change of variables Wj = uj — Uj-\ and Sj = Vj — Vj-i with 
the convention uq = vq = 0. For k = 1,2, ... ,m, define 



2m 2m 



Using the second inequality in (2.1), 

2m 2m 



CP 1 ">>>> -j £*ttlJ 

/ / $n,2m & y) exp { - -Var ( V ■ B%) - -Var ( V £ • ) dv d£ dy 

„ „ 2m 2m 

< C2 U |/fo)| II M* + «T P 
xEfnlexpf, -t %£L 



2m 

dw ds dy, 

- J=l 1 N W^j+i + s j+i n " + ' 

where ^/kJjXj (1 < j < 2m) are independent copies of the d-dimensional standard normal random 
vector and X2 m +i = 0. The rest of proof is similar to that of Proposition 3.3 in [3]. ■ 

Recall the definition of D in (3.10). For I = 1,2, ... , m and K > 0, define 

Dr/ = D Pi {u2i — U2e-i > K/n or V21 — v^t-i > K/n}. (3.13) 

The following result implies that the domain D 1 ^ £ does not contribute to the limit of even moments. 

Lemma 3.6 For any a £ @* and £ = 1, 2, . . . , m, 

r r 2 2m 

lim lim sup / / $ n , 2m (£, y) exp { - -Var ( TT & " <' f ' )) j d« dv d£dy = 0. 
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Proof. Let t = max(ti,t 2 ). Define 



^K,t = {u,V £ [0, t] 2m : Ui < U 2 < ■ ■ ■ < U 2m , Vl<V 2 <-'-< V 2 m, 

un ~ u 2 e~i > K/n, v 2 £ - V21-1 > K/n} 



and 

2m 



I n K,i = J R4md $n,2m(Z, v) e*P { - ^Var (n^i- - <'', )) } du dv di dy. 



'K,l j=l 

Applying Cauchy-Schwartz ineuqality, 



jn < 



(r r ..2m ^ i 

/ / *n,2m(e, y) exp { - -Var ( ]T ^ • (i^- 1 - B J 2 ) ) } du dv d£ dy 



According to Lemma 3.5, we can replace D^g in the above inequality with 

D K,l = D K,£ n | 2 U ^ ~ U2l ~ X ~ ' n ' 2 ^ ~ U2£_1 '"J ' 

The rest of the proof is similar to that of Proposition 3.4 in [3]. ■ 

We next divide into two subsets. In section 3.4, we will show that permutations in one subset 
do not contribute to the convergence of even moments. 

For each a in we introduce the following decomposition of/={l,2,..., 2m}: 

Jg e = {j G I : j is even and a(j) is even}, Ig = {j G I : j is even and o~(j) is odd}, 
loe = {j ^ I '■ J i s odd and a(j) is even}, J^ = {j £ I : j is odd and <r(j) is odd}. 

Let ig = {2j : 1 < j < m] and 7 = {2j - 1 : 1 < j < m}. Then I e = l% e + 1% and /„ = + l a 00 
for all (j in & . We make the change of variables w 2 k = u 2k , w 2 k-i = n{u 2k — u 2 k-i), s 2k = v 2k , 
S2fe-i = n(v 2 k - v 2 k-i) for k = 1, 2, ... ,m. Define 

L> n = |u>, s G M 2m : < w 2 < ^4 < • • • < w 2m < ii, < s 2 < s 4 < • • • < s 2m < t 2 

< w 2k -i < n(w 2k - w 2k -i), < s 2k -i < n(s 2k - s 2k -i), 1 < k < mj. (3.14) 

From the above decomposition of /, 

/2m 

= E { / n /(""K 1 - n w - - ""i^;, ..^ - 



x 



(3.15) 
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Assume that X2, X4, . . . , x 2m and Z2, 24, . . . , z^ m are linearly independent elements in some linear 
space. For any a in let 

^ee = { x j ~ z a(j) '■ 3 £ ^ee}> ^oe = { x j+l ~ z a(j) '■ 3 £ -foe}i 
^eo = { x j ~ z cr(j)+l '■ 3 £ Ieo}i = {Xj + i — Z CT (j)+l : 3 £ -Co}- 

Note that elements in each of the above sets are linearly independent. For simplicity, we use 
-j^A to denote the cardinality of a set A. Suppose = r - Then #A° e = #A^ = r and 

#A° = #A a oe = m-r. We are interested in the dimension of the set A a := A% e U A a oe U A£ U A£ D , 
that is, the maximum number of elements in A a which are linearly independent. Since elements 
in A a ee U A a oe are linearly independent, the dimension of A a is greater than or equal to m. 

Lemma 3.7 The dimension of A a is m if and only if {(J,<r(j)) ■ j € I% e } = {(j + l,<r(j) + 1) : 
j e r oo \ and {(j + l,a(j)) : j G I£J = (0>(j) + 1) : j € / e ff }. 

Proof. It suffices to show the only if part. Note that the m elements in ^4g e U are linearly 
independent. If one of the two condition fails, then there must exist an element in A a eo U A° Q such 
that it does not belong to the space spanned by A a ee U A a oe . This implies that the dimension of A a 
is greater than m. ■ 

Let £Pq = {o~ £ & : the dimension of A a is m} and = & — &q. Lemma 3.7 implies 

3.4 Convergence of even moments 

We will show the convergence of all even moments. Recall that 

/2m 
\{f{n H B^-n H B^ U) )dudv , 

J> j = l 

where 

D = jo = U < Ui < U 2 < ■ ■ ■ < U 2m < h, = V < Vi < V 2 < ■ ■ ■ < V 2m < ^j- 

Note that we can find a sequence of functions /at, which are infinitely differentiable with compact 
support, such that J Rd /at (a;) dx = and 



lim 

iV-KX) 



\f(x)-f N (x)\ (\x\^- d Vl)dx = 0. 



Therefore, by Proposition 3.1, we can assume that / is infinitely differentiable with compact 
support and J Rd f(x)dx = 0. 

We first show that permutations in V\ do not contribute to the limit of even moments using 
Lemmas 3.4 and 3.6. 

Proposition 3.8 For any a G 

ton^E [ / f{n H B% x ~ n H B^) dudv] = 0. (3.16) 

3=1 
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Proof. For any e > and K > 0, define 

Dk }€ = Dn{u 2 e-u 2 e-i < K/n, v 2 e-v 2 e-i < K/n, u 2i -u 2 i- 2 > e, v 2i -v 2 e- 2 > e, £=1,2,... ,m}. 

Thanks to Lemmas 3.4 and 3.6, we can replace D in (3.16) with D 1 ^ . 

Recall the equality in (3.15). Making proper change of variables gives 
2m 



11 



--E 



lJD K,ej=l 

{ L n / w - B ztfi n /w - - ^** 2 

"D^XT^ AtrTa Ar-J<J 



jel; 



3 el, 
x 



cr 
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H,2 



1) 



H\t>H,2 



S <r(j) + 1- 



dw; (is 



where 



o2m 



!»,s£ K"" : e < w 2 < W4 < ■ ■ ■ < w 2m <ti, e < s 2 < s 4 < • • • < s 2m < t 2 , 

< W2k-i < K A n(w 2 fc - w 2k _ 2 ), < s 2fc _i < K A n(s 2fc - S2/C-2), 

w 2k - w 2k _ 2 > e, s 2k - s 2 fc- 2 > e, fc = 2, ...,m|. (3.17) 

Since a G 'Pi, there exists a jo G 7f U i^, such that jo ^ l% e U J„ e . Without loss of generality, we 
can assume that jo G I° . Let p n (x,y) be the density function of the random field Z n = (Xj^J 1 ) 
where 



X; 



B^-B?f {j) if jG 

.<i-s5' if je/ CT e 



(3.18) 



and 



s 



Mj'o) + l ~ 



fiS-i] if J = Jo 



if j€l£. 



Since / is infinitely differentiable and had compact support, 

2m 



n 



m(2+Hd) 



E 



' e 5=1 



H\T}H,2 



f[n H [B«' L -B«>* }-n tl [B 1 
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H,2 



dw ds 



,mHd 



c 2 n 



c 2 n 



< c 3 n 



-Hd 



-Hd 



D'i 



: m +2+#ig e )d 



\f(n H x jo +y jo )\ \f(n H Xj)\ \f(n H x j +y j )\p n (x,y)dxdydwds 



[ m +2+#Ig e )d 



Pn{ u 1 V ) dx dy dw ds 



X 



D". JWL(l + #lSe)<i X 
K.e 



supp n (— 77, y) dy dw ds. 



n 
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Let Q n (w,s) be the covariance matrix function of Z n = (Xj,YJ l ) defined above. Then Q n has 
the following expression 

A C n 
C n B n 



Q, 



where A = A(w,s) is the covariance matrix function of X = {Xj), C n = C n (w,s) the covariance 
matrix function of X = {Xj) and Y = (Yj), and B n = B n (w, s) the covariance matrix function of 
Y = (Y j ). 

After doing some algebra, we have 



Qn 



(A — C^B n ^C n ) 1 —A l Cn(B n — C n A 1 C^') 1 

-B n C n (A — C n B n C n ) (B n — C n A C n ) 



and det(Q n ) = det(^4) det(I? n — C n A l C^). For simplicity of notation, we write 



D 1 Dl 
D 2 D A 



Note that 



x,y)Q n 1 (x,y) T = xDtx 1 + xD 2 y A + yD 2 x 1 + yD 4 y 1 
= xDix T + 2xD 2 y T + yD±y T 

= {x^/Dl^x^/Dif + 2(x^W l )(^/W l )- l Dly T + yD A y T 
> y(D 4 - D 2 D^Dl)y T . 



TT 



Then 



supp n (— w ,y) dy 



supp n (x,y) dy 



= 04/ (det(Q„)) 2 supexpj - \{x,y)Q n 1 {x,y) T \ dy 

< c 4 I (det(Q n ))-3 exp { - \y(D A - D 2 D^D^)y T \ dy 

= c 5 (det(Q„))-^(det( J D 4 - D 2 D^Dj))^ 
= c 5 (det(A))-2. 



Therefore, 



n 



m(2+Hd) 



E 



2m 



I) 



Hf(n H B^-n H B»f 3) )dudv 



< C6 n 



-Hd 



K,e j=l 



(det (A(w, s))) idwds. 



From the definition of X in (3.18), we see that the components of X are linearly independent 
and thus A(w, s) is not singular. Taking into account the definition of in (3.17) and the 

continuity of det(A(w,s)), we obtain det(A(w, s)) > c e > for all (w,s) in D\ e . Therefore, 



n 



m(2+Hd) 



E 



2m 



D 



Wf^B^-^B^dudv 



K,t j = l 



< CjC e 2 n 



This completes the proof. 

We next show the convergence of all even moments. 
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Proposition 3.9 For any m G N, 



lim E[F n (t 1 ,t 2 )] 2m = D% >d \\m d ^Wa{^t u t 2 )(\ 2m - 
Proof. By Proposition 3.8 and equation (3.11), we only need to show 



„ 2m 

Bmptajin^ e E [/„n /w -" hb S) 

o-G^n .7=1 



(in dt> 



!_ d E[Va(0,ti,t 2 )C] 



2i/f 



The proof will be done in several steps. 

Step 1 Since <r G by Lemma 3.7, the equation (3.15) can be written as 



n 2m E 



2m 



I) 



Hfi^B^-^B^dudv 



j'=i 



(3.19) 



E 



[ / ii /(-'«■' - <y) ii /(»'[«»■' - - »"^' rM - b&j) 



'(J)" 



We introduce random fields X n (w, s) = {X?(w, s) : j G J e } and y n (u>, s) = {Y^to, s) : j G I e } 
with 

if j€i£ 



^ R^' 1 _ R H < 2 



rR^' 1 _ R H > 2 1 



Mj')+i » 



B% +1 ], if i G 4- 



and 



if jG J e ° 



^[B 11 ' 1 ^ -B^]+n H [B 



H,h i Hfr>H,2 



Mj) -"S«r(j) + 1 



Let Z n {w,s) = (X n (w,s),Y n (w,s)). Denote the covariance matrix and the probability density 
function of the Gaussian random field Z n (w, s) by Q n (w, s) and 

p n (x,y) = (27r)- md (detQ n Ks))"5 e xp{ - -(x, y)Q n (w, s)" 1 ^, y) T }, 
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respectively. Then 



2m 



D 



3=1 



n mHd E 



n mHd 



On 



I f{n H X?(w, s))f(n H Xj(w, s) + Y?(w, a)) dw ds] 



j'6/e 



imd J JJ 



[ f(n H Xj)f(n H Xj + yj-i)p n (x,y)dw dsdxdy 



F(x, y) p n (—rr, y) dw ds dx dy, 
n n 



(3.20) 



where F(x,y) = Y\ f{x j )f(x j + yj). 

We need to compute the limit of the density p n [x/n H , y) as n tends to infinity. The covariance 
matrix between the components of X n (w,s) and Y n (w,s) converges to the zero matrix, and the 
covariance matrix of the random field Y n (w,s) converges to a diagonal matrix with entries equal 
to WjHi + when j G I° e and Wj^i + s afj) when j G I° Q . Let A cr (w,s) be the covariance 

matrix of X(w, s) = (Xj(w, s) : j G J e ) with 

Xj(w, s) = < 



ifiG/-. 

We see that the covariance matrix of the random field X n (w, s) converges to A a (w, s). Thus, 



lim p n (-g,y) = (2Tr)- md {detA°(w,s)) 



n + -sg)_i)~' «p ( - ^ 



^ ^-1 + S a(j)-1 



n 



( W 2H +S 2H 



1 l%f 



o ,.,2-ff , 2ff 
Z + S *(j) 



On the other hand, the region D n converges, as n tends to infinity, to 

jio, s G M+ m : < w 2 < w 4 < ■ ■ ■ < w 2m < h, < s 2 < s 4 < • • • < s 2m < 
Note that we can add a term —1 because f Rd F(x,y) dyj = for all j in I e and 



o Jo 



3 , i ivj-r 



oo />oo 



J 



'l _ g 2 W 2H +S 2H \ d w d s _ 



Therefore, provided that we can interchange the limit and the integrals in the expression (3.20), 
we obtain that the limit equals 



D 



HA 



d (27r) _2 ^ / (det A a (w,s)) ' 2 dwds, 



(3.21) 



o 
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where 

O = |0 < W 2 < U>4 < • • • < W2m < h, < S 2 < S 4 < • • • < S 2m < i 2 j. 

Finally, the left hand side of (3.19) equals 

D m r 
(M ! E -^r\\f\\t- d (^r md (det A* (w,s))->dwds 

= (2m-l)\\D% d \\f\\1 [ (27r)-^(detA(u,v))- 1 2d U dv, 

H a J E m 

and, taking into account of Lemma 2.2, this would finish the proof. 
Step 2 Recall the notation D n in (3.14). Define 

D n ,K = D n n |o < w 2 k-i < K A n(w 2k - w 2 k-2), < s 2 k-i < K A re(s 2fc - s 2 fe_ 2 ), 1 < fc < m|. 

The region D n ^ is uniformly bounded in n and we can then interchange the limit and the integral 
with respect to w and s, provided that we have a uniform integrability condition. 

Observe that 

/ \Pn(Ar>y)\ P dwds 

< ci / (det Qn{w, s)) 2 dsds 

j (! exp{-ivar(E^-^Ks) + E^-^K^))}^^ dwds 



C 2 



and 

Var S ) + E ^ • Y/>, s)) = Ii(£, t?) + J 2 (£, t?), 

where 

/iCe, 17) = Var ( E & • < A + E W • ^[^fi-fci " B ^) 

and 

Applying Cauchy- Schwartz inequality gives 

/ exp{-i[Ji(e,T;)+l2(^»/)]}ded»7 
<(/ expj-A^r?)}^^)^ / exp{- I 2 (£,7,)} 
Using similar arguments as in the proof of Proposition 3.4 in [3], we obtain 

„ m 

/ exp{ - ii (£,??)} d£dr] < c 3 ]7(u 2fc „i)~ Ha! (u 2A . ^— ^ - u 2fc „ 2 ) 



18 



and 



„ m 

/ exp { - I 2 (C,V)} di dr] < c 4 TT {v 2 k 
JR 2md t i 



\-Hdl v 2k-l \-Hd 

fe=i 



Therefore, for all p such that 1 < p < 



sup/ \pn(^jT,y)\ P dwds < c 5 , 



where C5 is a positive constant independent of n, x and y. 
Let 

Jn,if = / / F(x,y)p n (—n-,y)dwdsdxdy. 

Thus, taking into account that the function F(x,y) is continuous and has compact support, by the 
dominated convergence theorem, we obtain 

lim/ n ^= / F(x,y)( lim / Pn(Ar , y) dw ds) dx dy . 

On the other hand, there exists p > 1 such that 

f x 

SUp/ \Pn{— u-,y)\ P dwds < OO, 
" JD n>K n 

which implies 

lim I n K = / / F(x,y) lim 1 D ((w, s)) p n {—fj,y) dw ds dx dy. 
With the same notation as in Step 1 we get 

lim I n>K = / [det A°(w,s)]-hwds 



n— >oo 



X 



(2*) 

/ F(x,y) T / / (u 2H + v 2H y2(e~ 2 ^ m ^ m -l)dudvdxdy. 
Jm. 2md ~-r Jo Jo 



o 

K rK A |y,l 2 



3 eh 

The right hand side of the above equality converges to the term in (3.21) as K tends to infinity. 
Step 3 We need to show that 

lim limsup / / F(x,y) p n {—f7 ,y) dw ds dx dy = 0. (3.22) 

K^oo n^oo y R 2md J D n —D n j { 71 

Recall the equation (3.20) and the notation e in (3.13). 



x 

F(x,y)p n (^,y)dwdsdxdy 
smd -I D n —D n> K n 

2m 

n 



m(2+Hd) E r f J] f{n H B^ - n H B*f. ) du dv 



• /U « = l^if,£ J = l 

Therefore, the statement in (3.22) follows from Lemma 3.6. The proof is completed. ■ 

Proof of Theorem 1.1. This follows from Propositions 3.1, 3.3 and 3.9 by the method of 
moments. 
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4 Appendix 



Here we give some lemmas which are necessary in the proof of Theorem 1.1. 



Lemma 4.1 Assume that 1 < Hd < 2. There exists a positive contant c such that 



JO 



(w 2H + s 2H )~% dwds<c(aA b) 



2-Hd 



Proof. Without loss of generality, we can assume that a < b. Making the change of variable 
v = s/w gives 



a rb ra 

I (w 2H + s 2H )~2 dwds = / / 
) Jo Jo Jo 

pa /*oo 

< / / w^ Hd 
Jo Jo 



l ~ Hd {l + v 2H )-2 dvdw 
(1 + v 2H y^ dvdw 



where c\ is a positive constant independent of b. 



Lemma 4.2 Assume that 2 — H < Hd < 2. Let X be a d- dimensional centered normal random 
vector with covariance matrix a 2 1. Then, for any n E N and y S R rf , there exists a positive constant 
c depending only on H and d such that 



/ / ( w 2 H +s 2Hyi M 

Jo Jo 



exp I i 



y-X 



n H \/w 2H + s 2H 



dwds < cn Hd - 2 \y\TJ- d 



Proof. It suffices to show the above inequality when y ^ 0. Making the change of variables 
u = \y\~vnw and v = \y\~~Hns gives 



( w 2H + s 2H )-tE 



exp ( i 



y-X 



o Jo 

n Hd - 2 \y\^- d f°°(u 2H + v 2H )~% E exp (, 

Jo Jo \y\Vu 2H + v 2H 



n H Vw 2H + s 2H 

d 



dw ds 
y-X 



dudv 



<n Hd - 2 \y\i~ d 



/"OO c oo 

J J (u 2H +v 2H )-i(2A(u 2H + v 2H ym\X\)dudv 



= C n Hd - 2 \y\^- d . 

The last equality follows from using polar coordinates and the assumption 2 — H < Hd < 2. 
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